WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




PCT 

INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification 6 : 

H04B 15/00, H04K 1/00, H04L 27/30 



Al 



(11) International Publication Number: WO 99/40698 

(43) International Publication Date: 12 August 1999 (12.08.99) 



(21) International Application Number: PCT/US98/24367 

(22) International Filing Date: 13 November 1998 (13.1 1.98) 



(30) Priority Data: 

09/020,970 



9 February 1998 (09.02.98) 



US 



(71) Applicant: MOTOROLA INC. [US/US]; 1303 East Algonquin 

Road, Schaumburg, IL 60196 (US). 

(72) Inventors: STEWART, Kenneth, A.; Apartment 209, 5501 

Carriageway Drive, Rolling Meadows, IL 60008 (US). 
QIAO, Yi; 2931 South Wells Street, Chicago, IL 60616 
(US). ' 

(74) Agents: HAAS, Kenneth, A. et al.; Motorola Inc., Intellectual 
Property Dept., 1303 East Algonquin Road, Schaumburg, IL 
60196 (US). 



(81) Designated States: BR, JP, KR, European patent (AT, BE, CH, 
CY, DE, DK, ES, Fl, FR, GB, GR, IE, IT, LU, MC, NL, 
PT, SE). 



Published 

With international search report. 



(54) Title: 



METHOD AND APPARATUS FOR JOINT DETECTION OF DATA IN A DIRECT SEQUENCE SPREAD SPECTRUM 
COMMUNICATIONS SYSTEM 



z 



(401) 



KAS 
YQ 



R£<HV^ 4 °S?GNATURE 



ESTIMATES 



CHANNEL 
ESTIMATOR 



SEQUENCES (K,kJ 



convolution 
processor 

KAS 
YQ 



c 



(4(B) 



0ATA SYMBOLS 
DETECTOR 



(32V 



SOFT DECISION 
Iff 0RUATI0N TO ERROR 
CORRECTING 0EC00ER 



tUDAMBLE SEQUENCE 
m 
(404) 



SIGNATURE SEQUENCES 
(406) 



(57) Abstract 

Joint detection of data signals occurs in a code division, multiple access (CDMA) communication system as follows: a digital 
signal processor (315) first extracts a midamble portion of transmitted signal vectors and generates an estimate of the channel response 
corresponding to each user-antenna pair using a channel estimator (401). In the preferred embodiment, each user's transmission within the 
communication system comprises either the type-1 or type-2 burst with varying midamble and guard period durations. Next, a convolution 
processor (402) forms the convolution of the user signature sequence with the channel impulse response estimate associated with each 
user and each of the antennas (308-310). A detector (403) utilizes the set of vectors and creates a set of sub-system matrices and solves 
the sub-system matrices to extract symbol information from an individual remote unit's transmission. Finally, the detector (403) outputs 
symbol information for subsequent infrastructure processing. 
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METHOD AND APPARATUS FOR JOINT DETECTION OF DATA IN A 



Direct Sequence Spread Spectrum Communications System 



Field of the Invention 



5 



The present invention relates generally to joint detection of data signals and, in 
particular, to joint detection of data symbols simultaneously transmitted over a shared 
communications channel and distinguished by user-specific signature or spreading 
sequences. 



Since the early 1 SOT's, code division, multiple access (CDMA) cellular 
communications systems have been designed based on direct-sequence spread- 

15 spectrum (DS-SS) principles. In such systems, multiple users simultaneously occupy 
the same radio frequency channel, separated only by user-specific spreading or 
signature sequences. Probably the best known contemporary example of this 
approach to cellular system design is that defined by Personal Station-Base Station 
Compatibility Requirements for 1.8 to 2.0 GHz CDMA Personal Communication 

20 Systems" (American National Standards Institute (ANSI) J-STD-008), or, broadly 
equivalently, the Telecommunications Industry Association Interim Standard 95 (TIA IS- 
95). 

Recently, proposals for so-called M 3-rd Generation" cellular communication 
systems have been made to the European Telecommunications Standards Institute 



10 



Background of the Invention 
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(ETSI) for adoption within the Universal Mobile Telecommunications System (UMTS) 
standardization process. Among the proposals are systems based on Direct-Sequence, 
Spread-Spectrum (DS-SS) technology. One such candidate system, which is 
particularly relevant to the present invention, is usually referred to in the UMTS 
5 community as the TD-CDMA system (for "Time-Division Code-Division Multiple 
Access"). That system makes use of a combination of time- and code-division 
techniques as a means of improving overall system capacity while purportedly retaining 
compatibility with 2 nd Generation systems, notably the Group Special Mobile (GSM) 
cellular communications system. 
10 As can be seen from FIG. 1, the TD-CDMA air interface incorporates both TDMA 

and CDMA elements. The TDMA component is provided by dividing each radio 
frequency channel into frames of 4.615 ms duration, with each frame further divided 
into time slots of approximately 577us in length. The CDMA element is permitted by the 
allocation to each user of a unique 16-ary Walsh-Hadamard orthogonal spreading code 
15 which is used to spread the Quadrature Phase Shift Keyed (QPSK) or Quadrature 
Amplitude Modulation (16-QAM) data symbol sequences comprising the useful part of 
each user's transmission. 

The mobile station (MS) or base station (BS) transmission within any timeslot is 
referred to a 'burst 1 . As presently envisaged, the TD-CDMA air interface supports two 
20 distinct burst types, the general structure of which is shown in FIG. 1 . The so-called 
type-1 burst transmits a data sub-burst of 28 data symbols, followed by a midamble of 
length 296 chips (used for channel estimation purposes), a second sub-burst of 28 data 
symbols, and finally a guard period of 58 chips; the type-2 burst transmits 34 symbols 
in each sub-burst, with midamble and guard period durations of 107 and 55 chips 
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respectively. In both cases the length of each burst is 1250 chips. The same burst 
structures are used for both forward and reverse links, although since the forward link 
represents a point-multipoint transmission, the content of the midamble segment is 
slightly different in that case. This is of no consequence, however, for the present 

5 purposes. Accordingly, the description of the drawings below will focus on reverse link 
operation, pointing out differences between the forward and reverse only when these 
are significant for the description of the invention. 

FIG. 2 shows a communication system employing a TD-CDMA air interface. A 
single timeslot is shown by the figure, within which K mobile stations 201-209 are 

10 simultaneously active. In the example shown, the mobile station population 201-209 
transmits simultaneously on a specific timeslot, distinguished by length-16 spreading 
codes c (,) through c (K) respectively. It is important to note that - unlike IS-95 systems 
in which the sequence used to spread each symbol is a sub-sequence of a much 
longer sequence - the same code c ik) is used continuously to spread each data 

15 symbol from the same user. In other words, the received signal at BS 200 comprises a 
plurality of time overlapping coded signals from individual mobile stations, each 
^ transmitted within the same timeslot and distinguishable by a specific signature 
sequence. 

The use of conventional DS-SS receivers such as the RAKE receiver is not, 
20 however, envisaged for use in the TD-CDMA system. Rather, receiver designs capable 
of simultaneously or jointly recovering the data symbols transmitted by the population 
of mobile stations operating within the same timeslot/frequency within a given cell or 
sector are intended. Examples of such receivers are described in articles Lupas R., 
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Verdu S., "Linear Multiuser Detectors for Synchronous Code-Division Multiple-Access 
Channels", IEEE Trans. Inf. Theory, vol. 35, no. 1, Jan. 1989, Klein A., Baier P.W., 
"Linear Unbiased Data Estimation in Mobile Radio Systems Applying CDMA", IEEE J. 
Sel. Areas Comm., vol. 11, no. 7, Sept. 1993, Blanz J., Klein A., Nashan M., Steil A., 
5 "Performance of a Cellular Hybrid C/TDMA Mobile Radio System Applying Joint 
Detection and Coherent Receiver Antenna Diversity", IEEE J, Sel. Areas Comm., 
vol.12. No. 4, May 1994, and Jung P., "Joint Detection with Coherent Receiver Antenna 
Diversity in CDMA Mobile Radio Systems", IEEE Trans. Veh. Tech., vol. 44, no. 1, Feb. 
1995. 

10 The receivers described in the above-mentioned references perform joint 

detection using block linear joint sequence detectors. The computational complexity of 
the methods used to solve the system of equations arising from the above methods 
can, however, be very large. This results in receivers that are costly and consume large 
amounts of power. Therefore, a need exists for a method and apparatus for joint 

is detection of data signals that requires less computational complexity, resulting in less 
costly, lower power receivers. 

Put more generally, a need exists for a method and apparatus to reduce the 
computational complexity of linear joint detectors where the signature sequences of the 
users are constant over a block (or burst) of symbols, and where some external means 

20 are provided - in the form, for example, of a midamble or pilot sequence - for 
performing channel estimation. 
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Brief Description of the Drawings 

FIG. 1 illustrates a prior-art Time-division, Code-division Multiple Access (TD- 
CDMA) air interface incorporating both Time Division Multiple Access (TDMA) and 
Code Division Multiple Access (CDMA) elements. 
5 FIG. 2 is a block diagram of a communication system utilizing the TD-CDMA air 

interface of FIG. 1 for a single timeslot. 

FIG. 3 is a block diagram of a mobile station and base station illustrating 
transmission and reception of signals in the TD-CDMA communication system of FIG. 2 
in accordance with the preferred embodiment of the present invention. 
io FIG. 4 is a block diagram of the digital signal processor of FIG. 3 in accordance 

with the preferred embodiment of the present invention. 

Detailed Description of the Drawings 

Although the following description relates to the candidate UMTS TD-CDMA 
is system, it will be readily appreciated that the invention may apply to any CDMA system 
in which the signature sequences assigned to each user do not change from symbol to 
symbol over some fixed interval of time, and where some means of channel estimation 
are provided. For the particular TD-CDMA example, the interval is the burst (or more 
precisely, the sub-burst) and the midamble provides the means of channel estimation. 
20 Joint detection of data signals according to the invention in such a CDMA 

communication system proceeds generally as follows. A digital signal processor (DSP) 
first extracts the midamble portion of the burst (FIG. 1 ) received on each antenna and 
generates, using the channel estimation function, an estimate of the complex-valued 
channel impulse response defined from each user to each antenna. Next, the DSP 
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forms the convolution of the assigned user signature sequences with the channel 
impulse response estimate associated with each user and each antenna. A detector, 
as described by the detailed description of the invention below, then utilizes the 
resulting set of convolved signal vectors to create a system of sub-matrices and then 
5 solves that system of sub-matrices in order to extract the underlying modulated data 
symbol information from each user. The detector outputs soft-decision symbol 
information for use in subsequent error control decoding. Unlike prior-art detectors, 
creating the system of sub-matrices to extract the data symbol soft decision information 
results in less computational complexity than the direct approaches of the references 

10 mentioned above, resulting in receivers which are less costly and consume lower 
amounts of power compared to the prior-art joint detection receivers. 

Now, in more detail, Fig. 3 illustrates a general signal model within which the 
present invention may be described. In the preferred embodiment of the present 
invention each user {an example of which is mobile station 301) utilizes a particular 16- 

15 ary spreading sequence 303 to spread each of the symbols 302 transmitted within a 
burst (i.e. each data sub-burst of FIG. 1), although any spreading sequence length may 
be used. This is distinct from existing DS-SS cellular communications systems such as 
J-STD-008 or IS-95, in which both short and long codes have repetition intervals that 
are much greater than the symbol duration. Note that FIG. 3 shows only transmission of 

20 the first sub-burst of the burst of FIG. 1 since this is all that is necessary for description 
of the preferred embodiment. Transmission of the user midamble sequences and data 
sub-burst 2 is implicit. 

Returning to FIG. 3, during a specific timeslot, the first mobile station 301 
transmits a sequence of N complex-valued symbols 302 with each symbol spread 
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using spreading sequence c (,) 303. The spread signal is then filtered by a chip pulse- 
shaping filter p{t) 330, frequency-converted to carrier frequency /, using frequency 
converter 305, and transmitted by antenna 306. All the ^-l other mobile stations 
behave similarly, using codes c <2 > through 6 K) , respectively. 
5 BS 200 supports spatially diverse signal reception by providing K 0 widely 

separated antennas 308-310, which may be omni-directional or sectored antennas. 
Since BS 200 provides K a antennas, a set of KK a baseband-equivalent radio- 
frequency channels may be defined, with channel h (ik -\t,T) describing the complex 
channel impulse response linking mobile station k to base station antenna k a . 
10 (Assuming mobile stations have only a single antenna, the case of forward link 
transmission is obtained by simply setting K„ =i ). 

The signal received on base station antenna 308 is frequency converted to 
complex baseband form by mixer 311 and then processed by chip pulse matched filter 
314 before being sampled by a quadrature analog-digital converter (ADC) 317. The 
15 sampled signal is subsequently distributed to Digital Signal Processor (DSP) 320 
where the present invention is located. The signals received from other antennas 309- 
310 are similarly converted to sampled baseband signals for distribution to DSP 320. 
In the preferred embodiment of the present invention ADC 317 operates at the chip 
rate, i.e. at \/T c sample rate, although it will be obvious that the invention may be 
20 easily extended to include oversampling designs in which the ADC samplerate is 
pi T c where p is an integer oversampling rate. 
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In the preferred embodiment of the present invention mobile station transmit 
timing correction (similar to that employed in the GSM cellular communication system) 
is employed for the reverse link of the TD-CDMA system such that the bursts of the K 
simultaneously active users are observed quasi-synchronously at the BS receiver. 
5 That is, the bursts are received at the BS receiver aligned in time except for a timing 
error of the order of a fraction of the symbol duration T S1 or equivalent^ of the order of 

a few chips duration. 

Accordingly, ADC's 317-319 may be configured to sample (at chip spacing) the 
nominal received burst interval corresponding to the burst structure of FIG. 1, with any 
10 residual timing error incorporated simply as a shift in the channel estimate. As a result, 
data symbol recovery for each timeslot may be performed by examining the length- 
1250 (as stated above the timeslot duration is 1250 chips) complex-valued baseband 
signals s (kJ recovered from each antenna, and which span the duration of the received 

timeslot. The extent of the received signal s ( * -) is shown in FIG. 1 . 

15 FIG. 4 is a block diagram of DSP 320 in accordance with the preferred 

embodiment of the present invention. The blocks indicate specific functions to be 
implemented, mapped onto a programmable DSP with possibly multiple ALU's, or onto 
an application-specific VLSI device. In the preferred embodiment of the present 
invention DSP 320 would be a DSP such as a Texas Instruments TMSC320C80 or 

20 TMSC320C60 processor, or some other suitable programmable DSP. 

DSP 320 first extracts the midamble portion of the received signal vectors s ( * J 
and generates an estimate within channel estimator 401 of the channel h (k ' ka) (t y r) 
corresponding to each user-antenna pair by making use of the known midamble 
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sequence m (see FIG. 1)- This operation may be done using a variety of known 
channel estimation methods, including matched filtering, or periodic inverse filtering 
such as the approach described in Steiner B., Jung P., "Optimum and Sub-optimum 
Channel Estimation for the Uplink of CDMA Mobile Radio Systems with Joint 
5 Detection", by European Trans. Telecom., vol. 5, Jan.-Feb. 1994. The output of the 
channel estimator 401 is a set of KX a channel impulse response estimates, available at 

7^-spacing, of length- W and represented by the vectors: 

T 

h -tfr, ./r 2 ,~.h w ) W 

where T is the transposition operator. 
10 Next, in the convolution block 402, the DSP 320 forms the convolution of each 

known user signature sequence c a> with the channel impulse response estimate 
associated with each user-antenna pair. That is, the convolution block 402 forms the 
convolution of the channel estimate h (k > ka) with the each code vector c {k) (for 

k a = l 7 ... 9 K g ) to generate a set of KX a vectors: 



15 



T 



This set of KX Q vectors then enter the data symbol detector 403. 

In order to describe the invention, consider demodulation of only the length- 
A^data sequence comprising the first sub-burst of the burst structure illustrated in FIG 
1. Demodulation of the second sub-burst is performed in similar fashion and so need 
20 not be described separately. 
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Say the desired length-// transmitted data sequence for sub-burst 1 of user k 
is denoted by the vector rf A) : 



where the complex data symbols are taken from an M-ary alphabet V = {v P v 2 , -^v^) . 
5 Any alphabet may be used in theory, although in practice the TD-CDMA system 
proposes 4-ary (QPSK) or rectangular 16-ary (16-QAM) alphabets. 
Introducing a transmitted symbol vector as 



where 



10 By constructing a sub-system equation (will be described in detail later) as 



u„ 



i4fi,id„ + w„ 



(6) 



decoder 403 utilizes the sub-system to solve for d„ in order to recover the data symbols 
</<*\ w = l,...,Ar transmitted by MS k. 

The description of a preferred embodiment of the present invention will be 
15 outlined in three parts, covering: 

1 . Generation of the Sub-System Equation (Eq.(6)); 

2. Solving the Sub-System Equation for d„ Using a Forward Propagation 
Multiuser Joint Detection Technique (FPMJD); and 
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3. Solving the Sub-System Equation for d„ Using a Forward-Backward 
Propagation Multi-User Joint Detection (FBPMJD) technique. 

Generation of the Sub-System Equation 

5 

As discussed above, each of the data symbols w = l,...,tf transmitted by 
MS k is spread by user-specific signature sequence: 

c to={cp,cp.--.,c$>) i: (7) 

consisting of 0, generally complex-valued chips at chip interval T e . 
10 At antenna k a% the T c -spaced signal observation over the duration of the 

symbol sequence tf* } is denoted by vector e ( * fl) , of length nq+w-\. 

<*») (fa) (fa) (*> \ ka =l,-,K a (8) 

With a total of K a antennas, an overall, concatenated signal observation can be 
denoted by the vector e , 

T 

where e ( H ^=1";^ » s 9 iven in Eq.(8). 



15 In addition to the spread signal components, vector e also contains internally 
generated receiver thermal noise as well as co-channel interference due to users in 
different cells. These signal components are denoted by the vector n , 



T 

n 0) (2)" ... (K fl )-\ 



= ( 0) T <2) T ... (K«)Y 0°) 



where 
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Now, as a starting point in efficiently partitioning the detection problem, decoder 403 
defines a matrix partition parameter as, 

~Q + W-\ 



P = 



(12) 



s Decoder 403 then creates a set of sub-system matrices as: 



(to) 



(Ika) 



ilka) 
b m Q 



(2,to> 

b 


,(^.to) 


<m-l)0+l 




<2,to) 


(KM) 


h {m-\)Q+2 




(2,to) 


(KM) 











\<m<p (13) 



The p-th submatrix ^ is created by defining a zero-padding operator Ta defined as: 



<(2 - p)Q + W - 1) x <(2 - p)Q + W - 1) 



0 



Qx«2-p)Q+W-l) 



(14) 



and then defining the submatrix ^ created by decoder 403 as: 



(to)_ 

A P 1a 



<l,to) 


(2, to) 


(KM) 


(p-\)Q+i 


*(P-1X?+I 


b (p-\)Q+\ 


(Ika) 


(2,to) 


(KM) 






b (p-\)Q+2 


(Ilea) 




(KM) 









(15) 
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Similarily, decoder 403 creates a series of submatrices related to the last symbol 
vector j determination beginning from: 



and then using operator 



s to form 



(Ika) Qka) 
(Ilea) 9 

bQ< 



+ff'-l u Q+fV-\ 



(KM) 



(KM) 



(KM) 



Q*Q 



\q + w - 1) x (Q + w - 1) 



(16) 



(17) 



<*a> <*-> , m = 2,-,p 

AB,m lBAB,m-V 



(18) 



The sub-matrices ^,m=\-,p have dimension QxK , and j£,m=\—,p, have 

dimension (Q+W-\)*K, respectively. 

By constructing a Ka(NQ¥W-\)xKN matrix A as: 

^=(X>^> T >-->W") T ' < 19) 
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where 



A] 


0 


0 


0 


0 


0 


0 


K Ka* 
Al 


A\ 


0 


0 


0 


0 


0 


* ■* j 




A\ 


0 


0 


0 


0 




A p -\ 


A p -2 


ft 


\j -■ • 


A 
V 


n 
u 


0 


Ap 


A ko) ■ 

A p -\ 


0 


0 


o 


0 


0 


0 


}k Q > 

A P 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


A\ 


0 


0 


0 


0 


0 


0 


<*«> 

Ai 


}ka> 

A x - 


0 


0 


0 


0 


0 


}ka> 

A p -\ 




0 


0 


0 


0 


0 


' A P 


>„> 




0 


0 


0 


• 0 


0 


... 

A _ 


j4 _ _ 





(20) 



the received noisy signal vector e can then be presented by 
e=Ad+n (21) 
5 where d is given in Eq.(4), and e is given in Eq.(9), respectively. 

The system matrix A can be manipulated by the decoder 403 to make determination of 
the data symbol vector d more efficient, and so decoder 403 creates a set of matrices 
and vectors defined as 



(i) 



(2) 



<* fl > 



A B ., 



(1) 
AB,m 

(2) 
AB,m 



<*„> 

Ab^ 



(22) 
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where p , the matrix partition parameter, is given in Eq.(12), £f and are given in 
Eq.(13) through Eq.(18). 

In order to properly reflect the re-arrangement of the concatenated system matrix, the 
5 received signal vector e and noise vector n must correspondingly be re-arranged as: 

(23) 

r ( / a) 4'#«V + > - ««j?W-i] T 

(25) 

w ( / a) =l*w. - "feW-iP 

w„ = U I)T wfi» T - w^) T ] T . "V-" (26) 



allowing the modified data vector tn to be put in matrix form as 



10 
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HA m d n - m+] +A\d„ + w n 

m=2 
P 

m=2 
P 

S AB,m dN-m+l + AbJ d a/ + W n 



(27) 



Additionally by defining Un as: 



m=2 

r„~Z/lfl.mdw- m+ i " = A/ 



(28) 



the sub-system equation representing the relation between the transmitted symbols 
5 and the observed signals can be finally expressed as: 



n = l,~,N-\ 



Us.id n + w n • " mN 



(29) 



where A, and Ae , given in Eq.(22) form the corresponding sub-system matrices. 
In the embodiment of the present invention, two iteration procedures are proposed for 
estimating the transmitted data symbol vector d based on the sub-system form of 
Eq.(29); one is called forward propagation multiuser joint detection (FPMJD), the other 
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is called forward-backward propagation multiuser joint detection (FBPMJD). In the 
following sections, the two procedures are described separately in detail. 



Forward Propagation Multiuser 
Joint Detection (FPMJD) Technique 



Prior to describing the preferred embodiment of FPMJD, the following discussion 
sets out a necessary background for FPMJD analysis. 

Assume that jj m is a forward propagation estimate operator. When a zero- 
10 forcing (ZF) criteria is considered, M m is defined as: 



A#S = U H BA A B ^Ah (30) 



where ^ is used for estimating dn ,n = i,-,tf-i. while M ^ is used for estimating dN , 
respectively, 

According to Eq(29), if ^, is known, 4, can be obtained by 

M ™ u„ » = 1,-^-1 



A 



* ~" ' ' (31) 

Mb.zf u " n = N 



15 

However, as stated in Eq.(29), the vectors ^ /!=!,•■;#, are dependent on prior data 
symbols. By iteratively substituting the prior estimated data symbols 
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" ,m=n-m\n{n-],p-\},";n-\ t an estimate *,n=\— t N, of ^ can be obtained by 
decoder 403, i.e., 



A 



m-2 On-m+l 



m-2 
P 



a 

A 

OA/-J1 



B,m 

m ^2 UjV-m+1 



» = /> + l,.»,W-l (32) 
/i = A' 



The estimated ZF solution of the system in Eq.(29) can be obtained by 



(FP) A 

Mzf „ 

(FP) A 
M B.ZF 



/I = l,»-,tf-l 

n = N 



(33) 



10 



In the preferred embodiment of the present invention it can also be assumed that 
the noise and the transmitted symbols are uncorrected, i.e., the noise covariance 
matrix, denoted as R,,, and the transmitted symbol covariance matrix, denoted as 

are in the form of 

R n = a;I R d = a]l (34) 
where I is the identity matrix, cr and a) are the noise and symbol variances, 
respectively. 

When a minimum mean-square (MMSE) criteria is considered, the MMSE estimate 



operator jJ FF> is defined as: 
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2 " l 



U (FF) , (A?A>+?±1) A? 



d 



Similarity to the previously decided zero-forcing operator, A^^is used for estimating 

d ,« = i. -.N-i, and A^ /5e is used for estimating dN , respectively. The MMSE solutions 
of the system in Eq.(29) can be obtained by 



u " (36) 

1-.MMSE {FP) A n = ]sf 

M B.K1MSE „ 



In the preferred embodiment therefore, decoder 403 utilizes a FPMJD iterative 
procedure to estimate d,, as follows: 

A 

1) for n- 1 . . . jV, get * from r „ and estimated ^n^n-mir^n-Xp-l},--^^ (Eq.(32)) 

2) calculate the current estimation of ^ through * =a/ R) A (Eq.(33) or Eq.(36)) 

10 where M (FP) is the estimate operator defined in Eq.(30) or Eq.(35). 

This above described technique results in an estimate of * = with a slight 

loss of system performance due to the fact that detection of each symbol is based on a 
local and not global least-square or minimum-mean square criterion. To improve the 
system performance, in another embodiment of the present invention a technique 
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called forward-backward propagation (FBP) is utilized by decoder 403. This technique 
is described below. 

Forward-Backward Propagation 
Multi-User Joint Detection (FBPMJD) 

Prior to describing the preferred embodiment of FBPMJD, the following discussion 
sets out a necessary background for FBPMJD analysis. First , constructing, 



B= 



Ai 




Ab,i 






A2 


Ab,2 



(37) 



10 where, fafofa, and ^ 2 are given in Eq.(22) , and let 



W„-l 

w, 

/\ A 



2v, un Ad„ + A 2 d tv -i 



(38) 



(39) 



Zn 



, n = 2 t - 1 N 



(40) 



where w „,n=V" ^ is given in Eq.(25) and Eq.(26), and ^n=\~\K is given in Eq.(32), 
respectively. The vector ^ can be put in matrix form 



a„ 



n = N 



(41) 
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The system in Eq.(41) characterizes a backward propagation procedure for data 
symbol estimation, so the data symbol estimation operator of the system is called a 
backward propagation estimate operator denoted as }i^ BP > . 

Now if a zero-forcing (ZF) criteria is considered for data estimation, the 

5 backward propagation estimate operator K^ BP > can be defined as 

M (£P) = (B H B )~V 

where is used for estimating <k»=l-."-2, and is used for estimating d*_i, 

respectively, 

According to Eq(41), the backward propagation zero-forcing estimate of ^,n=v 
10 can be obtained by 



CBF) 

A 
A 

d n-l. ZF 



M zf ' (43) 

{BP) 



M B.ZF " 



Similarity, assume that the noise covariance matrix r„, and the transmitted 
symbol covariance matrix R<, are given in Eq.(34), then the backward MMSE estimate 
operator can be defined as 



lBP) (B"B +^4 1) B" 



era 

m7L e - WBbMvbV (44) 



15 
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and the backward propagation MMSE estimate of <k,n=i #-1, are obtained by decoder 
403 by solving: 



(BF) 

A 

A = 
&n-\ t MMSE 



Getting the solution from Eq.(43) or Eq.(45) is based on knowing vectors a n , which, in 
5 the preferred embodiment, can be obtained by decoder 403 through the forward 
propagation procedure (FP). That is, solving Eq.(43) or Eq.(45) actually involves two 
procedures: the forward propagation procedure characterized by Eq.(29) and the 
backward propagation procedure characterized by Eq.(41). We hence call the resulted 
solution of Eq.(41) the forward-backward solution. The detailed description to get the 

10 forward-backward solution a ,»=i,-,a/; through the forward-backward propagation 

procedure is presented as follows: 
For n=1,N 

A 

• get ^ from r„ and m=n-min {n-\ /?-!}• -*,/7-l . (Eq.(32)) 

• get a forward propagation estimation of (Eq.(33) or Eq.(36)) 

A 

15 • get z,, from ^ and estimated ^ m =„-\,n . (Eq.(39)) 

• get a,, from * and ^ . (Eq.(40)) 

Un~\ 



• get a backward propagation estimation of (Eq.(43) or Eq.(45)) 
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A 

A A 

• update the forward propagation solution: j = 

• End 

A 

• determine the last estimated data symbol: A = A 

Unlike prior-art estimation procedures, the above-described joint estimation techniques 

5 (FPMJD and FBPMJD) result in less computational complexity. The results in receivers 
which are less costly and consume lower amounts of power when compared to prior-art 
joint detection receivers. 

The descriptions of the invention, the specific details, and the drawings 
mentioned above, are not meant to limit the scope of the present invention. For 

10 example, although the above description was in the context of a reverse link, the 
invention may be applied equally to both forward and reverse links. Additionally, the 
above discussion referred to a generic burst structure comprising two distinct blocks of 
data symbols. Clearly both type-1 and type-2 TD-CDMA burst types are within the 
scope of the above discussion. Additionally, in order to support higher data rate 

15 services, it has been proposed that a single user may transmit on more than one 
timeslot with a frame, or may transmit more than one data symbol sequences (requiring 
more than a single spreading code to be allocated to that user). Nevertheless, it is 
intended that the present invention encompass any data symbol sequence that may be 
utilized. It is the intent of the inventors that various modifications can be made to the 

20 present invention without varying from the spirit and scope of the invention, and it is 
intended that all such modifications come within the scope of the following claims. 
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What is claimed is: 

1. An apparatus for detection of data signals from a transmitted signal in a 
communication system, the transmitted signal comprising a plurality of time overlapping 
coded signals transmitted from individual remote units, each coded signal 
distinguishable only by a specific signature sequence, the apparatus comprising : 

a channel estimator having a known signal as an input and outputting a channel 
estimate of a channel existing between a remote unit and a receiver antenna within the 
communication system; 

a convolution processor having the channel estimate and the signature 
sequence of the remote unit as an input and outputting a set of vectors representative 
of the transmitted signature sequence modified by a channel; and 

a detector having the set of vectors as an input and outputting symbol 
information for the remote unit's transmission, the symbol information being a solution 
to a set of sub-system matrices derived from a decomposition of a system matrix. 

2. The apparatus of claim 1 wherein the set of sub-system matrices are represented by 
the equation 



wherein u„ is the vector comprising the received signals and some prior estimated data symbols, 
Ai and A B ,i are subsystem matrices , w„ is noise, and d„ is the data symbol to be estimated. 

3. In a communication system utilizing Direct-Sequence, Spread-Spectrum (DS-SS) 
techniques for the transmission of signals, a method for determining an individual 
user's transmitted data symbols from a received signal, the method comprising the 
steps of: 




representing the received signal as a system matrix; 
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partitioning the system matrix into sub-matrices to form a plurality of sub-system 
matrices; 

solving the plurality of sub-system matrices via an iterative process to determine 
individual user's data symbols from the transmitted signal; and 

transmitting the individual user's data symbols to infrastructure equipment for 
further processing. 

4. The method of claim 3 wherein the step of solving the plurality of sub-system 
matrices further comprises the step of solving the sub-system matrices via an iterative 
process where an estimation of a transmitted symbol is based on an estimation of prior 
transmitted symbols. 

5. The method of claim 4 wherein the step of solving the plurality of sub-system 
matrices further comprises the step of solving the sub-system matrices via an iterative 
process where an estimation of a transmitted symbol is additionally based on 
subsequently transmitted symbols. 



BNSDOCID; <WO. 



.9940698A1J_> 



WO 99/40698 



PCT/US98/24367 



-26- 

6. In a Time-Division, Code-Division Multiple Access (TD-CDMA) communication 
system, where a transmitted signal comprises a plurality of time overlapping coded 
signals transmitted from individual remote units, each coded signal transmitted within a 
same timeslot and distinguishable only by a specific encoding, a method for 
determining an individual user's transmitted data from a received signal, the method 
comprising the steps of: 

extracting a midamble portion of the transmitted signal; 

convolutionally processing the midamble portion of the transmitted signal to 
produce a set of vectors representative of a received signature sequence; 

representing the convolutionally processed midamble portion of the transmitted 
signal as a system matrix; 

partitioning the system matrix into sub-matrices forming a plurality of sub-system 

matrices; 

solving the plurality of sub-system matrices via an iterative process to determine 
the individual user's data symbols from the transmitted signal; and 

transmitting the individual user's data symbols to infrastructure equipment for 
further processing. 
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